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A Possible Spin- Orbit Misalignment 
in the Transiting Eccentric Planet HD 17156b* 



00 

O 
O 
(N 

X) 



6 



(N 
> 

in 
o 



Norio Narita,^** Bun'ei Sato,^ Osamu Ohshima,^ and Joshua N. Winn" 
^ Department of Physics, School of Science, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan 

narita@utap.phys. s. u-tokyo. ac.jp 
^ Global Edge Institute, Tokyo Institute of Technology, 2-12-1 Ookayama, Meguro, Tokyo 152-8550, Japan 
^ Mizushima Technical High School, 1230 Nishiachi, Kurashiki, Okayama 710-0807, Japan 
* Department of Physics, and Kavli Institute for Astrophysics and Space Research, 
Massachusetts Institute of Technology, Cambridge, MA 02139, USA 

(Received 2007 December 16; accepted 2008 January 30) 
Abstract 

We present simultaneous photometric and spectroscopic observations of HD 17156b spanning a transit 
on UT 2007 November 12. This system is of special interest because of its 21-day period (unusually long 
for a transiting planet) and its high orbital eccentricity of 0.67. By modeling the Rossiter-McLaughlin 
effect, we find the angle between the sky projections of the orbital axis and the stellar rotation axis to be 
62° ±25°. Such a large spin-orbit misalignment, as well as the large eccentricity, could be explained as the 
relic of a previous gravitational interaction with other planets. 
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1. Introduction 

More than 250 exoplanets are now known, and their 
orbital characteristics are remarkably diverse. In partic- 
ular, the surprising discoveries of close-in giant planets 
and planets on highly eccentric orbits have led to inter- 
esting revisions of planet formation theory. It is generally 
believed that the close-in giant planets originally formed 
at larger orbital distances (beyond the "snow line") and 
migrated inward during the planet formation epoch (Lin 
et al. 1996). One possible migration mechanism is tidal 
interaction with the protoplanetary gas disk, after the 
planet opens up a gap in the disk (Type II migration; 
see, e.g., Lin & Papaloizou 1985), which is appealing be- 
cause it can account for the observed period distribution 
(see, e.g., Ida & Lin 2004). However, some planets have 
large eccentricities despite their close-in orbits, such as 
HAT-P-2b (Bakos et al. 2007), GJ 436b (Butler et al. 
2004), and HD 17156b (Fischer et al. 2007). Since disk- 
planet interaction would not excite a planet's eccentric- 
ity to such a level (Sari & Goldreich 2004), the scenario 
does not provide an obvious explanation for these systems. 
On the other hand, some other migration models would 
naturally produce eccentric orbits, namely, planet-planet 
gravitational scattering (Rasio & Ford 1996; Marzari & 
Weidenschilhng 2002; Nagasawa et al. 2008). In these 
models, planets would initially obtain very large eccentric- 
ities and small periastron distances through the gravita- 
tional interactions, and subsequently evolve into shorter- 
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period and more circular orbits through tidal dissipation. 

Assuming that planets are initially formed on circular 
orbits, the observation of highly eccentric orbits is already 
an indication of previous gravitational interactions. What 
other type of evidence might establish the case for such 
interactions? Recently, the alignment angle of the stel- 
lar spin axis and the planetary orbital axis (the spin-orbit 
alignment angle) has become recognized to be a promising 
diagnostic. This is because disk-planet interaction would 
probably maintain the original spin-orbit alignment, while 
planet-planet scattering would cause significant misalign- 
ments in a nontrivial fraction of cases (e.g., Chatterjee 
et al. 2007; Nagasawa et al. 2008). Because the damping 
time scale of the spin-orbit alignment angle is expected 
to be longer (by a few orders of magnitude) than that of 
eccentricity (Hut 1981; Winn et al. 2006), one might re- 
gard the spin-orbit angle as a fossil record of planetary 
migration. 

It is possible to measure the sky-projection of the spin- 
orbit alignment angle. A, for a transiting exoplanet, by 
making use of the Rossiter-McLaughlin effect (hereafter 
the RM effect: Rossiter 1924, McLaughlin 1924). The RM 
effect is the radial velocity (RV) anomaly caused by the 
partial occultation of the rotating stellar surface (see Ohta 
et al. 2005, Gaudi & Winn 2007 for theoretical descrip- 
tions), and several observers have reported the detections 
of the RM effect and measured A in transiting planetary 
systems (Queloz et al. 2000; Winn et al. 2005; Winn et al. 
2006; Narita et al. 2007; Winn et al. 2007; Loeillet et al. 
2007; Wolf et al. 2007). However, the reported angles for 
the spin-orbit alignment were either small or consistent 
with zero, even for the highly eccentric case of HAT-P-2b 
(Winn et al. 2007; Loeillet et al. 2007). 
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In this Letter, we present a photometric and spectro- 
scopic study of a recent transit of HD 17156b (Barbieri 
et al. 2007). This planet has a large orbital eccentricity 
and a small periastron distance, and is therefore a promis- 
ing candidate for a spin-orbit misalignment. The system 
has also attracted considerable interest because it has by 
far the longest period (21 days) of any known transiting 
exoplanet. We describe our observations in Section 2 and 
our results in Section 3. The final section is a summary 
and a discussion of the results. 

2. Observations and Data Reduction 

We observed a transit of HD 17156b on UT 2007 
November 12, both spectroscopically and photometri- 
cally, with telescopes in the Okayama prefecture of Japan. 
The transit was predicted to occur in the later part 
of that night, according to ephemerides kindly pro- 
vided by D. Charbonneau (2007, private communica- 
tion). The photometric transit was observed with a 
20 cm Meade telescope at Kurashiki-shi. The spectro- 
scopic transit was observed with the 188 cm telescope at 
the Okayama Astrophysical Observatory (OAO). We ob- 
served HD 17156 for a 5 hr period spanning the transit 
time, with both telescopes, through air masses ranging 
from 1.3 to 1.6. 

2.1. Photometry 

The photometric observations were conducted at 
Kurashiki-shi (E133°40'15", N34°32'27", h=7m) using 
a 20 cm Meade LX200R-20 f/10 telescope. A cooled 
(253 K) ST-9XE CCD camera provides a 17.6 x 17.6 
arcmin^ field of view with the pixel scale of 2.06 arcscc. 
We defocused the stellar images in order to average over 
pixel-to-pixel sensitivity variations, and to draw out the 
exposure time and thereby increase the duty cycle. The 
exposure time was 60 s and the readout time was 1 s. 
We observed through a Rc band filter covering the range 
5720-6780 A. 

We used the Astronomical Image Processing for 
Windows software (AIP4Win Ver.2.1.10) for the subtrac- 
tion of dark current, flat-fielding, and aperture photome- 
try. We determined the apparent magnitudes of HD 17156 
(y = 8.17) and two comparison stars: HD 16906 (F = 8.28) 
and BD+71 168 {V = 9.57) using an aperture radius of 20 
pixels. The typical FWHM of each star ranged from 18 
to 21 pixels (from 37 to 43 arcsec). We estimated the 
sky background level with an annulus from 27 to 32 pixels 
in radius centered on each star, and subtracted the esti- 
mated sky contribution from the aperture flux. Then we 
computed the differential magnitude between HD 17156 
and the ensemble average of the comparison stars. After 
these steps, we clearly detected the transit event. We also 
found small linear trend in the out-of-transit (OOT) data; 
we removed the trend by fitting the OOT data with a lin- 
ear function of time. 

We estimated photometric errors of our samples as fol- 
lows. We first calculated the Poisson noise of each sample 
and found that the scatter for the OOT data is systemat- 
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Fig. 1. Top: Photometric liglit curve observed witli the 
20 cm telescope (251 samples) on UT 2007 November 12. 
The error-bars are scaled to satisfy X^/'^dof = l-O- Bottom: 
Radial velocity samples on that night computed from the 
OAO/HIDES spectra. The values and uncertainties are pre- 
sented in Table 1. 

ically larger than the Poisson error. In order to account 
for the systematic errors in the photometry, we scaled the 
Poisson estimates of the errors to satisfy x^/'^dof = 1-0 
for the OOT dataset. Next, we also considered the time- 
correlated noise (the so-called "red noise", see e.g., Pont 
et al. 2006). We compared the standard deviation of each 
OOT sample (cri) with the standard deviation of the OOT 
data after averaging each N points (ctat). Wc found that 
aN < (Ti/\/N for values of N ranging from 10 to 60, sug- 
gesting that the effects of correlated noise are small in 
this case. Thus we did not modify the error bars further. 
The typical photometric error is about 4 mmag. Our light 
curve is shown in the top panel of Figure 1. 

2.2. Spectroscopy and Radial Velocity 

We used the High Dispersion Echelle Spectrograph 
(HIDES; Izumiura 1999) on the 188 cm telescope at 
OAO. We employed the standard setup for RV measure- 
ments with the Iodine absorption cell (covering the wave- 
length range 5000-6200 A, where the iodine absorption 
lines are abundant). The slit width of 280 /im (1'.'06) 
yielded a spectral resolution of ^45000, and the see- 
ing was between 1'.'5 and 2'.'0. For the transit event of 
UT 2007 November 12, we used 15 min exposures with a 
binnningof4 (spatial) x 1 (dispersion). We also obtained 
a few spectra on UT 2007 November 14, 17, and 18, out- 
side of transits, in order to refine the spectroscopic orbit. 
For these OOT spectra we used 30 min exposures and 1 
X 1 binning. The typical signal-to-noise ratio (SNR) was 
50-80 pixel^^. We processed the frames with standard 
IRAF^ procedures and extracted one-dimensional spec- 
tra. We then calculated relative RV variations by the 
algorithm following Sato et al. (2002). We estimated in- 
ternal errors of the radial velocities from the scatter of 

^ The Image Reduction and Analysis Facility (IRAF) is dis- 
tributed by the U.S. National Optical Astronomy Observatories, 
which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Table 1. Radial velocities obtained with the OAO/HIDES. 



Table 2. Star and planet parameters. 



Time [HJD] Value [m s '-] Error [m s '-] 



2454417.11678 


94.98 


25.53 


2454417.13164 


108.60 


16.24 


2454417.14398 


119.28 


16.71 


2454417.16256 


78.87 


19.01 


2454417.17275 


99.93 


15.94 


2454417.18389 


75.77 


16.24 


2454417.19495 


90.34 


16.74 


2454417.20621 


64.90 


16.37 


2454417.21594 


59.26 


15.77 


2454417.22968 


84.55 


21.55 


2454417.23773 


45.40 


19.12 


2454417.25048 


54.91 


14.57 


2454417.26162 


39.53 


16.10 


2454417.27249 


6.80 


13.14 


2454417.28340 


18.77 


13.29 


2454417.29502 


21.32 


19.48 


2454417.30549 


5.71 


15.86 


2454417.31692 


-6.74 


14.82 


2454417.32796 


33.70 


14.45 


2454417.33911 


28.36 


15.27 


2454417.34993 


5.32 


16.25 


2454419.08778 


-140.28 


10.18 


2454421.92713 


36.65 


10.99 


2454422.17592 


66.31 


18.30 


2454423.13311 


88.78 


12.96 



the RV solutions for 4-5 A segments of the spectra. The 
typical errors are 10-20 m s~^. Our radial velocities are 
presented in Table 1, and plotted in the bottom panel of 
Figure 1. 

3. Analyses and Results 

As described in the previous section, we have obtained 
25 RV samples and 251 Rc band photometric samples 
covering the transit. In addition, in order to determine 
an optimal solution of orbital parameters for HD 17156, 
we incorporate previously published radial velocity data 
taken at the Subaru telescope and the Keck I telescope (9 
from Subaru and 24 from Keck: Fischer et al. 2007). 

We employ the analytic formulas for light curves and ra- 
dial velocities including the RM effect given in Ohta et al. 
(2005) and Ohta et al. (2006) (hereafter the OTS formu- 
lae) in order to model the observed data. We also adopt 
the correction for Vsinlg in the OTS formulas by modify- 
ing V^sin/3(0TS) =Fsin/,(real)*l.l; see Winn et al. 2005 
and Narita et al. 2007 for details. (We note that this cor- 
rection does not affect the results for A, since in this case 
A and Vsinls are nearly uncorrelated parameters.) We 
adopt the stellar mass Mg = 1.2 [Mq] and the stellar ra- 
dius Rs = 1.47 [Rq] (Fischer et al. 2007), and the quadratic 
limb-darkening coefficient iti = 0.31 and U2,r ~ 0.35 for 
the Rc band photometry, and g = 0.49 and U2,g = 0.28 
for the spectroscopic band, based on Claret (2004). The 
adopted parameters are summarized in Table 2. 



Parameter 


Value 


Source 


A/. [Mq] 


1.2^ 


Fischer et al. 2007 


Rs [Rq] 


1.47 


Fischer et al. 2007 


Ul,r 


0.31^ 


Claret 2004 




0.35 =* 


Claret 2004 




0.49 =* 


Claret 2004 


a [A.U.] 


0.28^ 


Claret 2004 


0.15^ 


Fischer et al. 2007 


P [days] 


21.2162 ±0.0036 


This work.'' 




3.13±0.21 


This work.*^ 


[-^Jup] 


1.21±0.12 


This work.*^ 



Determined in this work thanks to D. Charbonnoau 2007, 
private communication. 

'^:Determined from the values and errors in this letter and 
the error of Ma and Ra presented in Fischer et al. 2007. 



Our model has 11 free parameters in total. Seven pa- 
rameters for the HD 17156 system include the RV ampli- 
tude K, the eccentricity e, the longitude of periastron 
the sky-projected stellar rotational velocity Fsin/s, the 
sky-projected angle between the stellar spin and the plan- 
etary orbit axes A, the ratio of star-planet radii Rp/Rs, 
and the orbital inclination i. We also add three param- 
eters for velocity offsets to the respective RV dataset vi 
(for our template spectrum) , V2 (for the Subaru data) and 
V3 (for the Keck data) , and one parameter for the time of 
mid-transit Tc on UT 2007 November 12. 

As a first step, we determined Tc using our photomet- 
ric data only, and used this result to refine the estimate 
of the orbital period P, with reference to a previous tran- 
sit epoch (D. Charbonneau 2007, private communication). 
We found Tc = 2454417.2645 ± 0.0022 [HJD], and thereby 
P = 21.2162 ± 0.0036 [days]. We adopt P = 21.2162 for 
the subsequent analysis (the uncertainty is negligible for 
our purpose). 

Next, we determined the orbit of HD 17156 by simul- 
taneous fitting of the photometric data and the RV data. 
Our fitting statistic is 



E 



i=l 

sin/, -2.6 



-\ 2 



2 W/=251 _ -,2 

Jj,obs Jj,calc 



E 



0.5 



(1) 



where Vcaic and fcaic represent the values calculated by 
the OTS formulae with the above parameters. The last 
term is a priori constraint on Vsinls, which enforces the 
spectroscopic determination by Fischer et al. (2007). We 
found optimal orbital parameters by minimizing the 
statistic using the AMOEBA algorithm (Press et al. 1992), 
and estimated confidence levels by Ax^ as the parameters 
were stepped away from the optimal values. (We also 
estimated the parameter uncertainties with bootstrap and 
Markov Chain Monte Carlo methods, as described below.) 

To account for possible systematic errors in the RV 
measurements (from photospheric jitter or instrumental 
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Fig. 2. Radial velocities of HD 17156 as a function of orbital 
phase, and the best-fitting model with the a priori constraint 
on y sin 7s. Three symbols represent the diff'erent dataset 
respectively (circle: OAO, triangle: Subaru, square: Keck). 
Left panel: The entire orbit. Right panel: A zoom of transit 
phase. Bottom panels: Residuals from the best-fit curve. 

sources), we used the following procedure. First, we de- 
termined the optimal parameter set for the OOT data 
using the nominal RV errors, and calculated the from 
each RV dataset (OAO, Subaru, and Keck). We then com- 
puted x^/i'dofi where lydof is the number of each dataset 
minus 4, which is the number of related parameters for 
the RV data {K, e, uj, and each RV offset). Next, when 
x'^/'^dof > 1, we inflated the radial velocity errors to sat- 
isfy /vdof = 1.0. After these steps, we re-calculated 
optimal parameters and uncertainties using all datasets. 
We also computed the results for the x^ statistic without 
the a priori constraint on Vshxlg, to support our claim 
that the application of this constraint does not affect the 
results for A. The resultant parameters are presented in 
Table 3^, and the RV curve with the best-fit model is 
shown in Fig 2. Consequently, our result for the key pa- 
rameter A was A = 62° ± 25° (60° ± 21°) with the reduced 
X^ of 0.97 (0.96). The numbers in parentheses refer to 
the case without the a priori constraint. The value of A 
is fairly large, and is inconsistent with zero at the ^2.b<7 
level, indicating a possibly large spin-orbit misalignment. 

Since the statistical significance of the result is mod- 
est, we checked on the calculation of the confidence levels 
in two different ways: a bootstrap analysis, and a Markov 
Chain Monte Carlo algorithm. For the bootstrap analysis, 
we first computed the residuals from the best-fitting model 
for UT 2007 November 12 data. Then we scrambled the 
residuals with their error-bars in a random manner, and 
created a new RV dataset by adding the residuals back 
to the best-fitting velocities. In this way, we created 1000 
fake data sets, and we calculated the optimal parameters 
in each case, using the same method that was applied to 
the actual data. The mean values and standard devia- 
tions of the 1000 results were A = 67° ± 26°, which are in 
excellent agreement with the x^ analysis. For the Markov 
Chain Monte Carlo algorithm, we used a variant of the 
code that has been employed by Holman et al. (2006) and 
Winn et al. (2006) for the Transit Light Curve project 
(see those papers for details, or Tegmark et al. 2004 for 
an introduction to the method). The resulting a posteri- 

Note that we present only Vsin/s and A as the results without 
the a priori constraint in the table 



Table 3. Fitted values and uncertainties^ of the parameters. 



Parameter 


Value 


Uncertainty 


K [m s^^] 


274.7 


±3.5 


e 


0.675 


±0.005 


u; [°] 


122.0 


±0.6 


Rp/Rs 


0.0846 


±0.0026 


^[°] 


85.65 


±0.29 


f 1 [m s "'"] 


146.7 


±3.9 


V2 [m s"^] 


92.9 


±2.6 


n [m s-i] 


93.4 


±1.2 


Tc - 2454417 [HJD] 
Vsinls ^ [km s^"'^] 


0.2645 


±0.0021 


2.8 


±0.5 


AMI (Ax^) 

A [°] (bootstrap) 

A [°] (MCMC) 


62 
67 

65 


±25 
±26 
±25 


ysin/s [km s^^] 


4.7 


±1.5 


A - [°] 


60 


±21 



=':Computed by Ax^ = 1.00. 

'':With the a priori constraint on Vsinlg. 

'^:Without the constraint. 



ori probability distribution for A was nearly Gaussian in 
shape, with mean 65° and standard deviation 25°. 

4. Summary and Discussion 

In this Letter, we have presented the results of si- 
multaneous spectroscopy and photometry of a transit of 
HD 17156b, which was reported to have an eccentric or- 
bit. (The high eccentricity is indeed confirmed by our 
radial velocity data around and after the transit phase.) 
We have measured the sky-projected spin-orbit alignment 
angle A by modeling the Rossiter-McLaughlin effect, and 
found A ~ 62° ± 25°. Although the statistical significance 
is modest, this result suggests that HD 17156b has a large 
spin-orbit angle. 

One may wonder whether the values e = 0.67 and A = 62° 
are physically possible and consistent with theoretical pre- 
dictions. Considering that the periastron distance q for 
HD 17156b is q = a{l - e) ~ 0.05 AU, where a is the 
semimajor axis and e is the eccentricity, and assuming 
a planetary tidal quality factor of Qp ~ 10^, then the 
damping timescale for the eccentricity due to dynamic 
tides in the planet is typically longer than 10 [Gyr] (Rasio 
et al. 1996; Ivanov & Papaloizou 2004). Thus, tidal cir- 
cularization is expected to be incomplete at its stellar age 
(5.7t} g [Gyr]: Fischer et al. 2007). There is another pos- 
sibility that the planet initially obtained larger eccentric- 
ity and more close-in periastron distance (say, e > 0.99 
and q ~ 0.04 [A.U.]) and evolved to the current orbit 
(Nagasawa et al. 2008). In this case, the planet would 
be still in the evolutionary track toward a hot Jupiter. 

On the other hand, the damping of the spin-orbit an- 
gle is mainly caused by tidal dissipation inside the host 
star, rather than the planet. Assuming a stellar tidal 
quality factor of Qs ~ 10^ as a typical value, the damp- 
ing timescale of the spin-orbit alignment is estimated to 
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be longer than ~ 10^^ [yr] (see e.g., Eq. (15) in Winn 
et al. 2005). We note that the tidal quality factors for 
the star and the planet are still unknown and the damp- 
ing timcscalcs estimated above remain fairly uncertain. In 
particular, it is possible that the coplanarization timescale 
is shorter than we have calculated, if only a thin outer 
convective layer is reoriented. However, our simple con- 
siderations show that it is at least physically plausible for 
the planet to have maintained a large eccentricity and a 
large spin-orbit angle since the planet formation epoch. 
Moreover, very recently several theoretical simulations 
predicted that certain degree of planetary systems mi- 
grated through the planet-planet interactions would have 
a large eccentricity and also a large spin-orbit misalign- 
ment (Chatterjee et al. 2007; Nagasawa et al. 2008). Our 
results are thus consistent in principle with the predictions 
from these recent theoretical works. 

Finally, we note that this planet is the first plausi- 
ble candidate to have a large spin-orbit misalignment. 
Further RV measurements are highly desired to bolster 
the precision in the determination of A. Such observa- 
tions are perfectly feasible with a larger telescope than 
the 1.9 m telescope we employed; the main difficulty is 
the rarity of follow-up opportunities from a given obser- 
vatory, due to the long orbital period. Given the present 
results, we believe this system and its transits are worthy 
of further photometric and spectroscopic scrutiny. 

We are very grateful to Yuuki Moritani, Akira Imada, 
Daisaku Nogami, and members of Okayama Planet Search 
Project for kindly exchanging OAO/HIDES observing 
time. We also thank David Charbonneau for providing 
us the previous transit time, and Shigeru Ida and Yasushi 
Suto for helpful discussions. We appreciate the careful 
reading and useful comments by the referee. Fred Rasio. 
N.N. is supported by a Japan Society for Promotion 
of Science (JSPS) Fellowships for Research (DC2: 18- 
10690), and B.S. is supported by Grant-in- Aid for Young 
Scientists (B) No. 17740106. 
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